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Abstract 

Consider a multi-user OFDMA system where multiple users share multiple discrete subcarriers, 
but at most one user is allowed to transmit power on each subcarrier. To adapt fast traffic and channel 
fluctuations and improve the spectrum efficiency, the system should have the ability to dynamically 
aUocate subcarriers and power resources to users. Assuming perfect channel knowledge, two formu- 
lations for the joint subcarrier and power aUocation problem are considered in this paper: the first is 
to minimize the total transmission power subject to quahty of service constraints and the OFDMA 
constraint, and the second is to maximize some system utihty function (including the sum-rate utihty, 
the proportional fairness utihty, the harmonic mean utihty, and the min-rate utihty) subject to the total 
transmission power constraint per user and the OFDMA constraint. In spite of the existence of various 
heuristics approaches, httle is known about the computational complexity status of the above problem. 
This paper aims to fiU this gap, i.e., characterizing the complexity of the joint subcarrier and power 
aUocation problem for the multi-user OFDMA system. It is shown that both formulations of the joint 
subcarrier and power allocation problem are strongly NP-hard as the number of users increases. The 
proof is based on a polynomial time transformation from the so-called 3-dimensional matching problem. 
Several subclasses of the problem which can be solved to global optimality or e-global optimality in 
polynomial time are also identified. 
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I. Introduction 

With the rapid growth of wireless data traffic, multi-user interference has become a ma- 
jor bottleneck that limits the performance of the existing wireless communication services. 
Mathematically, we can model such an interference-limited communication system as a multi- 
user interference channel in which a number of linearly interfering transmitters simultaneously 
send private data to their respective receivers. Exploiting time/space/frequency/code diversity 
are effective approaches to mitigate/manage multi-user interference. For instance, when the 
transmitters and/or receivers are equipped with multiple antennas, a joint optimization of the 
physical layer transmit-receive beamforming vectors for all users can efficiently mitigate multi- 
user interference [fT|-[|3|; when all transmitters and receivers are equipped with a single antenna, 
one way to control/mitigate multi-user interference is to impose certain frequency restrictions or 
transmission power limits on the frequency resources used by each transmitter Q, [[5|. 

Conventional Orthogonal Frequency Division Multiple Access (OFDMA) schemes preassign 
subcarriers to users in a nonoverlapping way, thus users (transmitting on different subcarriers) 
cause no interference to each other. Although the OFDMA scheme is well suited to be used in 
a high-speed communication context where quality of service is a major concern, it can lead 
to inefficient bandwith utilization. This is because that the preassignment of subcarries can not 
adapt traffic load and channel fluctuations in space and time; i.e., a subcarrier preassigned to a 
user can not be released to other users even if it is unusable when the user's channel conditions 
are poor. To adapt these fluctuations and improve the overall system's throughput, OFDMA 
based subcarrier allocation networks such as Worldwide Interoperability for Microwave Access 
(WiMAX) Q and Long Term Evolution (LTE) [VJ should be equipped with dynamic subcarrier 
and power allocation algorithms. In particular, a dynamic OFDMA based subcarrier and power 
allocation algorithm is well suited for the dense femtocell downlink system [[8j, where a large 
number of femtocells close to each other are deployed in a macrocell. 

The joint optimization of subcarrier and power allocations for the OFDMA system is a 
nonconvex problem, therefore various heuristics approaches have been proposed for this problem 
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Q' 0' [[8l-p3|. Among them, the reference [12| proposed a centralized OFDMA reuse parti 



tioning algorithm. The authors in [|8| developed a self-organized approach for OFDMA cellular 
networks, where each cell independently tunes its subcarrier and transmission power allocation. 
Recently, the authors in [|5| proposed a dynamic joint frequency and transmission power allocation 
scheme called enhanced Dynamic Frequency Planning (eDFP). eDFP is a hybrid centralized and 
distributed architecture, where a central broker first dynamically partitions subcarriers among 
neighboring cells so that the long-term cell-edge inter-cell interference is minimized, and then 
each cell independently allocates subcarriers and transmission power to its users in a way that 
its total transmission power is minimized. 

Notice that nonconvex optimization problems are difficult to solve in general. However, not all 
nonconvex problems are hard to handle since the lack of convexity may be due to an inappropriate 
formulation, and many nonconvex optimization problems indeed admit a convex reformulation 
||4|, [14 1, [15|. Therefore, nonconvexity of the joint subcarrier and power allocation problem for 



the multi-user OFDMA system does not imply that it is computationally intractable (strongly 



NP-hard in terms of computational complexity theory [16|). This paper aims to characterize the 
computational complexity of the joint subcarrier and power allocation problem for the multi- 
user OFDMA system; i.e., to categorize when the problem is (strongly) NP-hard and when it is 
polynomial time solvable. 

In fact, the dynamic spectrum management problem without the OFDMA constraint has been 
extensively studied in Q, [|9|, p7|-[19|. It is shown in Q, ||9| that the dynamic spectrum 



management problem is (strongly) NP-hard when the number of subcarriers is greater than two, 
or when the number of users is greater than one. However, the analysis of these results is 
highly dependent on the crosstalk channel gains among users; i.e., some of the crosstalk channel 
gains are assumed to be large enough, and some of them are assumed to be zero. We shall see 
late (in Section [II]) that the magnitude of crosstalk channel gains has no influence on the user's 
transmission rate if all users are required to transmit power in an orthogonal manner. This makes 
the two problems (the dynamic spectrum management problem without the OFDMA constraint in 
Q and the joint subcarrier and power allocation problem with the OFDMA constraint considered 
in this paper) sharply different from each other. An interesting result shown in [[9| is that the 
optimal solution of the two-user sum-rate maximization problem is automatically OFDMA if the 
crosstalk channel gains of the two users on each subcarrier are large enough. In addition, based 
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on the fact that the duality gap of the dynamic spectrum management problem is zero as the 
number of subcarriers goes to infinity (regardless of the convexity of the objective function) 



[T7[, p8| , dual decomposition algorithms have been proposed in p8| , p9| for the nonconvex 
optimization problem of multi-carrier systems, again without considering the OFDMA constraint. 

In this paper, we focus on the characterization of the computational complexity status of the 
joint subcarrier and power allocation problem for the multi-user OFDMA system. In particular, 
we consider two formulations of the joint subcarrier and power allocation problem. The first 
one is the problem of minimizing the total transmission power in the system subject to all 
users' quality of service constraints, all users' power budget constraints per subcarrier, and the 
OFDMA constraint (at most one user is allowed to transmit power on each subcarrier). The 
second one is the problem of maximizing the system utility (including the sum-rate utility, the 
proportional fairness utility, the harmonic mean utility, and the min-rate utility) while the total 
transmission power constraint of each user, individual power constraints on each subcarrier, and 
the OFDMA constraint are respected. The main contribution of this paper is that we prove the 
previously mentioned two optimization problems are strongly NP-hard, and we also identify 
several subclasses of the joint allocation problem which can be solved to global optimality or 
e-global optimality in polynomial time. 

II. Problem Formulation 

In this section, we introduce the system model and problem formulation. Consider a multi-user 
multi-carrier OFDMA system, where there are K users (transmitter-receiver pairs) sharing 
subcarriers. Throughout the paper, we assume that N > K; i.e., the number of subcarriers is 
greater than or equal to the number of users. Otherwise, the OFDMA constraint is infeasible. 

Denote the set of users and the set of subcarriers by /C = {1, 2, K} and J\f = {1, 2, A^}, 
respectively. For any k E IC and n E Af, suppose G C to be the symbol that transmitter 
k wishes to send to receiver k on subcarrier n, then the received signal at receiver k on 
subcarrier n can be expressed by 

where h^-EC is the channel coefficient from transmitter j to receiver k on subcarrier n and 
G C is the additive white Gaussian noise (AWGN) with distribution CN'{Q,r]^). Denoting the 
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power of by p^; i.e., pi := , the received power at receiver k on subcarrier n is given by 



where al^^ := stands for the channel gain from transmitter j to receiver k on subcarrier 

n. Treating interference as noise, we can write the SINR of receiver k on subcarrier n as 

Vk + EjVfc »k,jP] 
and transmitter fc's achievable data rate Rk as 



SINR" = , keIC, neU, 

K __r) I \ ^ _ n ..T) ' ' " 1 



i?fc = J]log(l + SINR^), keK. (1) 

neAT 

In this paper, we consider the following joint subcarrier and power allocation problem for the 
multi-user multi-carrier OFDMA system 



s.t. Rk > 7fc, k e )C, (2) 

Pfc > Pfc > 0, keK, neAf, 

pIp] = 0, y J ^k, neAf, 
where the objective function means the total transmission power of all users on all subcarriers, 
7a,. > (A; G /C) is the desired transmission rate target of user k, and p^ is the transmission power 
budget of user k on subcarrier n. The last constraint in ([2]) is the so-called OFDMA constraint, 
which implies that at most one user is allowed to transmit power on each subcarrier. 
Due to the OFDMA constraint, we know for any k e K, and n E M 

'~ rk + E,^k^hPl~^^' 

In fact, if p1 = 0, the above equality holds trivially; while if p^ > 0, it follows from the OFDMA 
constraint p^p" = (Vj ^ k) that p" = (Vj 7^ k) and thus '^j^^'^kjP] = O5 which shows 
that the above equality holds as well. Thus, problem (|2]) is equivalent to 



S.t. log (1 + ^) > ^ e 

Pfc > Pfc > 0, keIC, neAf, 

pIp^ = 0, W J ^k, neAf. 
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By introducing a group of binary variables x1 (A; G /C, n G A/"), problem ([2]) can be reformu- 
lated as 



mm 
s.t. 



(4) 



Vk 

xlPl > Pfc > 0, A; G /C, n G AT, 
G {0, 1} , kelC, ne Af, 

where the binary variable = 1 if user k transmits power on subcarrier n, or = otherwise. 
The last constraint J^keJC < 1 G A/") stands for the OFDMA constraint. 

Problem (|4]) can be dealt with by a two-stage approach. At the first stage, we solve the 
subcarrier allocation problem; i.e., determining the binary variables {x^}, which is equivalent 
to partitioning the set of subcarriers Af = {1,2,..., N} into K nonoverlapping groups {Mk}k=i- 
At the second stage, we solve the power allocation problem; i.e., solving K decoupled power 
allocation problems 

min > p2 

l^fe/nGATj. neA/'fe 

Pk > K > 0, ne Uk. 
Problem (|5]) at the second stage is convex, and thus is easy to solve. 

To sum up, the joint subcarrier and power allocation problem can be equivalently formulated 
as (|2]), (|3]), or (|4]). Formulation Q is intuitive and is easy to understand, whereas formulation 
Q is compact and is easy to analyze. The analysis of this paper is mainly based on Q. 

Besides the total transmission power minimization problem, we also consider the utility 
maximization problem for the multi-user multi-carrier OFDMA system, which can be expressed 
by 

max H(Ri,R2,...,Rk) 

s.t. J^Pk^Pk^ ^ ^ ^' 

neM (6) 

Pfc > Pfc > 0, keIC, neAf, 
PkPl = 0, ^ J 7^k, neAf, 
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k=l 



where H{Ri, R2, ...,Rk) denotes the system utility function and pk {k G JC) is the power budget 

of transmitter k. Four popular system utility functions are 

1 ^' 

- Sum-rate utility: Hi{Ri, R2, Rk) = —/^Rk 

k=l 

- Proportional fairness utility: H2{Ri, R2, Rk) = I JJ^-Rfe 

/ ^ \ 

- Harmonic mean utility: H3{Ri, R2, Rk) = K/ ^-R^^ , 

\fc=i / 

- Min-rate utility: if4(i?i, i?2, -Ri^:) = min [Rk] ■ 
It is simple to see that 

Hi{Ri,R2, Rk) > H2{Ri, R2, Rk) > -f^sl-Ri, R2, Rk) > Hi{Ri,R2, Rk), 
and the equality holds true if and only if Ri = R2 = ... = Rk- 

III. Hard Cases 

In this section, we show that both the total power minimization problem ([3]) and the system 
utility maximization problem ([6]) are intrinsically intractable (strongly NP-hard in the sense of 
computational complexity theory [jl6[|), provided that the ratio of the number of subcarriers and 
the number of users, that is N/K, is equal to any given constant number c > 1. To begin with, 
we briefly introduce complexity theory in Subsection III-A Then we show problems ([3]) and ([6]) 



are strongly NP-hard in Subsections |III-B| and |III-C[ respectively. 



A. A Brief Introduction to Complexity Theory 

In computational complexity theory p6| , a problem is said to be NP-hard if it is at least as 
hard as any problem in the class NP (problems that are solvable in Nondeterministic Polynomial 
time). The NP class includes well known problems like the 3-colorability problem (which is to 
check whether the nodes of a given graph can be colored in three colors so that each pair of 
adjacent nodes are colored differently). NP-complete problems are the hardest problems in NP 
in the sense that if any NP-complete problem is solvable in polynomial time, then each problem 
in NP is solvable in polynomial time. The 3-colorability problem is NP-complete. 
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A problem is strongly NP-hard (strongly NP-complete) if it is NP-hard (NP-complete) and it 
can not be solved by a pseudo-polynomial time algorithm. An algorithm that solves a problem 
will be called pseudo-polynomial time algorithm if its time complexity function is bounded 
above by a polynomial function related to the length of the given data of the problem. The 
3-colorability problem is also strongly NP-complete. However, not all NP-hard (NP-complete) 
problems are strongly NP-hard (strongly NP-complete). For instance, the partition problem is 
NP-hard but not strongly NP-hard. 

Strongly NP-hard or NP-hard problems may not be in the class NP, but they are at least as 
hard as any NP-complete problem. It is widely believed that there can not exist a polynomial 
time algorithm to solve any NP-complete, NP-hard, or strongly NP-hard problem. Thus, once 
an optimization problem is shown to be NP-hard, we can no longer insist on having an efficient 
algorithm that can find its global optimum in polynomial time. Instead, we have to settle with less 
ambitious goals, such as finding high quality approximate solutions or locally optimal solutions 
of the problem in polynomial time. 

The standard way to prove an optimization problem is NP-hard is to establish the NP-hardness 
of its corresponding feasibility problem or decision problem. The latter is the problem to decide 
if the global minimum of the optimization problem is below a given threshold or not. The output 
of a decision problem is either true or false. The feasibility or decision version of an optimization 
problem is usually in the class NP. Clearly, the feasibility or decision version of an optimization 
problem is always easier than the optimization problem itself, since the latter further requires 
finding the global minimum value and the minimizer. Thus, if we show the feasibility or decision 
version of an optimization problem is NP-hard, then the original optimization problem must also 
be NP-hard. 

In complexity theory, to show a decision problem B is NP-hard, we usually follow three steps: 
1) choose a suitable known NP-complete decision problem A; 2) construct a polynomial time 
transformation from any instance of A to an instance of B; 3) prove under this transformation 
that any instance of problem A is true if and only if the constructed instance of problem B is 
true. Furthermore, if the chosen NP-complete problem A is strongly NP-complete, then problem 
B is strongly NP-hard. 
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B. Strong NP-Hardness of Power Minimization Problem (|3]) when N/ K = c > 1 

To analyze the computational complexity of problem ([3]), we consider its feasibility problem. 
If the feasibility problem is strongly NP-hard, so is the original optimization problem. 

Feasibility Problem of (|3]). Given a set of transmission rate levels 7^, individual power 
budgets per subcarrier p'^, direct-link channel gains a^f., and noise powers t]^, check whether 
there exists a subcarrier and power allocation strategy such that 

^ log 1 + > 7fc, e /C, 

' Pl>p'^>0, k e IC, neAf, ^'^^ 
= 0, V j ^ fc, neAf. 

To analyze the computational complexity of the feasibility problem (|7]), we choose the fol- 
lowing strongly NP-complete problem (see [16]). 

3 -Dimensional Matching Problem with Size K. Given three different sets X, y, and Z 
with |A:'| = |3^| = |Z| = _ft', and a relationship set TZ, where each triple corresponds to a 3-way 
marriage that would be acceptable to all three participants. The problem is to check whether 
there exists an arrangement of these K marriages so that polygamy is avoided and everyone 
receives an acceptable spouse. 

We shall show that any instance of the 3-dimensional matching problem corresponds to an 
instance of the transmission rate feasibility problem 

Theorem 3.1: Checking the feasibility problem of (|3]) is strongly NP-hard when N/K = 2. 
Thus, problem ([3]) itself is also strongly NP-hard when N/K = 2. 

Proof: Assume that N/K = 2. Consider any instance of the 3-dimensional matching 
problem with 



X = 


11 2 




y = 




■,Ky}, 


z = 







and a relationship set 

^ = {(kxjyjz) I k^ E x,jy ey,ke z}, 
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we construct a multi-user multi-carrier system where there are K users (which correspond to set 
X) and 2K subcarriers (which correspond to set y [j Z). More exactly, IC = X and M = y[J Z. 
Define 

^2 = {{h,lz) \ {K,jy,lz) e TZ} . 

For each k E /C, the power budgets per subcarrier p^, the noise powers 77^, and the direct-link 
channel gains f, are given by 



(8) 




1, 
2, 
3, 



and 



ni no 



if {k, n) G iSi; 
if (fc,n) G ^2; 
if (A;,n) ^5iU'52, 

1, if (fc,ni,n2) e 7^; 
0.25, if (A;,ni,n2) ^ 7^, 



respectively. Letting 



7fc = 3 log 2, A; e /C, 



the corresponding instance of problem Q is 



1 + 



> 3 log 2, ke /C, 



3 > K > 0, A; G /C, n G 3^, 
2 > K > 0, A; G /C, n G 
p^p^ = 0, M J ^k, neM, 



(9) 



(10) 



(11) 



(12) 



(13) 



where 77^ and a^j^. are given in ( fTO| ) and ( fTT[ ), respectively. We are going to show that there exists 
an arrangement such that everyone is happy (i.e., the answer to the 3-dimensional matching 



problem is yes) if and only if the constructed problem ( [13] ) is feasible. 

We first show that if the answer to the 3-dimensional matching problem is yes, then problem 



( fT3| ) is feasible. In fact, if {{kx,jy,lz)} is an appropriate arrangement of the marriage, a feasible 
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power allocation of problem ( [T3| ) is given by 



3, 
2, 
0, 



ifn = jy; 
if n = Iz] 
if n 7^ or h. 



(14) 



This is because, since {{k^, jyjz)} is an appropriate arrangement of the marriage, the above 
power allocation strategy ( [T4| ) is orthogonal to each other (i.e., Pl^p^^ =0, ^ jx, ^ n E M). 
Furthermore, we have for user k^ = lx,'^x, ■■■,Kx, 

f^x ^r^x^ f^x 



Rkx 



(from ([13])) 




(from ^ 



(15) 



(from (|9]) - ([TT]) and ([T4j)) 



3 log 2. 



So ( [T4] ) is a feasible power allocation of problem ( [T3] ). 



On the other hand, we show that if all the constraints in ( 13) are satisfied, the answer to the 3 



dimensional matching problem must be yes. Notice that for any user k^ E /C, if the transmission 
rate Rkx > 3 log 2, it must transmit on at least two subcarriers, for otherwise if it transmits only 
on one subcarrier, its maximum transmission rate is 

1*3' 

■ 2 1og2<3 1og2. 



log 1 + 



Since there are K users and = 2K subcarriers and at most one user is allowed to transmit 



on each subcarrier (by the OFDMA constraint), the feasibility of ( [T3] ) asks that each user in the 
network must transmit on exactly two subcarries. By the construction of the parameters of the 
system, we can verify that the corresponding direct-link channel gains of the user on the two 
subcarriers must be 1, for otherwise the transmission rate is at most 



log 1 



1*3\ , / 0.25*3 
+ log 1 + 



^ . , -o , - . ^ . log7< 3 log2. 

Furthermore, the fact that all users' transmission rate requirements are satisfied implies that each 
user transmits on one subcarrier in {ly,2y, Ky} with noise power 1 and one in {1^, 2^, i^^} 
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with noise power 2, which makes the transmission rate equal to 

, / 1*3\ , / 1*2\ 

log 1 + —— + log 1 + -— = 3 log2. 



Otherwise, 

- if one user transmits on two subcarries in {l^, 2^, K^} , the transmission rate is at most 

(1 * 2\ / 1 * 2\ 

1 + — j +log h+—j =2 1og2<3 1og2. 

- if one user transmits on two subcarries in {ly, 2y, Ky} both with noise power 1, then at 
least one user will transmit on two subcarries in {1^, 2^, K^} either with noise power 2 
or 3 and the transmission rate of this user is at most 2 log 2. 



Therefore, problem ( [13] ) is feasible if and only if that for all users kx {k^ G /C) there exist jy 
and L such that 



«fe' k =(^k k 



1, ikx,jy)&<Sl, {kx,lz)^S2 (16) 



and 

K 



{ly, 2„ Ky, 1„ 2„ K,} = U {jy, / J . (17) 



k=l 



According to the construction of Si, S2 and ( [TT| ), we see that 

is an arrangement of the happy marriage. 

It is simple to check that the above transformation from the 3-dimensional matching problem 



to the feasibility problem ( 13 ) can be done in polynomial time. Since the 3-dimensional matching 



problem is strongly NP-complete, we conclude that checking the feasibility of problem ( [T3| ) is 
strongly NP-hard. Therefore, the optimization problem ([3]) is also strongly NP-hard. ■ 
To illustrate the above proof, we take the following 3-dimensional matching problem as an 
example, 

{Ix; 2a;, 3a;, 4a:;} ; 

y = {Ij/, 2y, 3y, 4y} , 

^ — {^zy'^zy'^zy^z} : 
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and 



It is simple to check that 

h)} (18) 

is a perfect match for the above given instance of the 3 -dimensional matching problem, since 
every one receives an acceptable spouse. Based on this 3-dimensional matching problem, we con- 
struct a 4-user 8-carrier system with /C = {l^;, 2^, 3^, 4a.} and Af = {ly, 2y, 3^, Ay, 1^, 2^, 3^, A^}. 
According to ([8]), we have 

'^1 = {(l^', 2y), (2a;, ly), {S ^ , 2 y) , {S ^ , Ay), {A^, 3y)} 

and 

S2 = {{U, 2.), (1., 4,), (2,, 2,), (2,, 3,), (3,, 2,), (3a;, 3,), (4,, 1,)} . 

The proof of Theorem |3.1| suggests the following system parameters (cf. (|9l)-([TT])). All power 
budgets per subcarrier are 

Pk 3, k la:, 2a;, 3a;, 4a;, ^ 1?^' 2y) 3j^, 4y, 
Pk 2, k Ij., 2a., 3a;, 4a., TL 1^, 2-,, 3^, 4^, 

all noise powers are r^^ = 3 (/c G /C, n G Af) except 

2ii l-u 2-1/ 4^1 3i/ 1 

Vil = V2I = Vsl = Vsl = vZ = 1, 
ft'?'' = rif^ = 772. = = r7o" = ?7o" =77!" =2; 

'-La; ' ' ' ^x '^x '^x '^x ' 

and all direct-link channel gains are a^j^. = 0.25 {k E )C,n E M) except 



a 



a 



2y 

^X j^X 



a 



Ix 1 la; 



a 



at", 

^X 1-^X 



a 



2a: ^^x 



3a; i3a; 



a 



3a; )3j3 



3a: ;3a: 



2x )2x 

3y 

C^A A 
^X -.^x 



a 



2a; i2a; 



a 



1. 



In this example, based on the solution (18), we can construct an OFDM A solution to the 



corresponding feasibility check problem (13), i.e., 



PU 5 PU ' Plx ' Plx ' Plx ' Plx ' PU ' Pla 



P2x ' P2x 5 P2x ' P2x ; P2x ' ^'2^ 5 F2^ ) F2, 



Lj, 2„ 3y Ay I, 2, 3z 4, 

P-i , P-i , P3 , P3 , P-i , 5 , P3 

Ox ' ■'■ Ox ' ■'■ Ox ' ^ Ox ' ^ Ox ' ^ Ox ' ^ Ox ' ^ -Ja 



/'4^ 5 ^4:^ 5 /'4;, 5 F4:i; 5 F4^ 5 ^4:^ 5 P ix ' /^4, 



(0,3,0,0,0,0,0,2), 
(3,0,0,0,0,2,0,0), 
(0,0,0,3,0,0,2,0), 
(0,0,3,0,2,0,0,0). 



(19) 
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On the other hand, to look for a feasible solution of problem ( [13] ), we have to make each 
user transmit on two subcarriers (one with noise power 1 and the other with noise power 2) 
and the direct-link channel gains of the user on the corresponding two subcarriers be 1. Notice 
that r]ll = r]ll = r]ll = r]ll = 1, r]^: = vt = vt = vll = 2, and = a^^^^ = 

al^ n = al^ 9 = at" o = a,^ o = a!" 4 = a]'' , =1. We can ask user l^^ to transmit 

i^X I'^X Ox i^X Ox ^X i^X ^X i^X 



on subcarriers 2y and 4^, user 2^ to transmit on subcarriers \y and 2^, user 3^ to transmit on 
subcarriers Ay and S^, and user 4^ to transmit on subcarriers 3^ and 1^, respectively. Consequently, 
{(la;,2y,42), (2a;, ly,2^), (33;,4y,3^), (42:,3y, I2)} Is di solution to the 3-dimensional matching 
problem. 

Theorem |3.1| shows that checking the feasibility of problem Q is strongly NP-hard when 



c = 2. Based on this basic result, we can further prove that it is strongly NP-hard to check the 
feasibility of problem ([3]) when N /K = c provided that c is a strictly greater than one constant. 

Theorem 3.2: Given any constant c > 1, checking the feasibility of problem (|3]) is strongly 
NP-hard when N/K = c. Thus, problem ([3]) itself is also strongly NP-hard. 



Proof: By Theorem 3.1 we only need to consider the cases: (i) 1 < c < 2 and (ii) c > 2. 
In case (i), we partition K users into (c — 1)K Type-I users and (2 — c)K Type-II users, and 
= cK subcarriers into 2(c — 1)K Type-I subcarriers and (2 — c)K Type-II subcarriers. We 
construct the channel parameters of Type-I users on Type-I subcarriers in the same way as in the 
case that c = 2. Furthermore, noise powers, direct-link channel gains, and power budgets of all 
Type-II users on Type-II subcarriers are set to 0.3, 1, and 3, respectively; these parameters of all 
Type-I users on Type-II subcarriers and all Type-II users on Type-I subcarriers are set to 3, 0.25, 
and 1, respectively. All Type-II users' desired transmission rate targets are set to be log 11. See 
Fig. [T] for the corresponding multi-user multi-subcarrier OFDMA system. Our construction is 
such that the channel condition of the Type-I (Type-II) user on the Type-I (Type-II) subcarrier 
is reasonably better than the one of the Type-I (Type-II) user on the Type-II (Type-I) subcarrier. 
We can check that the only possible way for all K users to meet their transmission rate targets is 
that, each Type-II user transmits full power on one Type-II subcarrier (actually Type-II users and 
subcarriers are dummy users and subcarriers) and all Type-I users appropriately transmit power 



on Type-I subcarriers. By Theorem 3.1[ however, checking whether all Type-I users' transmission 



rate requirements are satisfied is strongly NP-hard. 
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Fig. 1. An illustration of the multi-user multi-carrier OFDMA system constructed for case (i) and case (ii), where the number 
of Type-I subcarriers is exactly twice of the number of Type-I users. In both cases, the channel parameters of Type-I users 
on Type-I subcarriers are chosen in the same way as in the case that c = 2. Noise powers, direct-link channel gains, power 
budgets of Type-II user(s) on Type-II subcarriers are set to 0.3, 1, and 3, respectively; these parameters of Type-I users on 
Type-II subcarriers and Type-II user(s) on Type-I subcarriers are set to 3, 0.25, and 1, respectively. For case (i), the transmission 
rate targets of Type-II users are all set to log 11. For case (ii), the transmission rate target of the single Type-II user is set to 
((c-2)Jf + 2)logll. 



In case (ii), we partition K users into K — 1 Type-I users and 1 Type-II user, and = cK 
subcarriers into 2{K — 1) Type-I subcarriers and {c — 2)K + 2 Type-II subcarriers. We construct 
the channel parameters of Type-I users on Type-I subcarriers in the same way as in the case 
that c = 2. Moreover, noise powers, direct-link channel gains, and power budgets of the single 
Type-II user on Type-II subcarriers are set to 0.3, 1, and 3, respectively; these parameters of 
all Type-I users on Type-II subcarriers and the single Type-II user on Type-I subcarriers are 
set to 3, 0.25, and 1, respectively. The transmission rate of the single Type-II user is required 
to be not less than ((c — 2)K + 2) log 11. See Fig. [T] for the corresponding system. Due to 
special construction of the system, we can check that the only way for all K users to meet 
their transmission rate targets is that, the single Type-II user transmits full power on all Type-II 
subcarriers and all Type-I users appropriately transmit power on Type-I subcarriers. Again, by 



Theorem 3.1 checking whether all Type-I users' transmission rate requirements are satisfied is 
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Strongly NP-hard. ■ 

C. Strong NP -Hardness of Utility Maximization Problem ([6]) when N/ K = c > 1 

This subsection considers the utility maximization problem (|6]) for the multi-user OFDMA 
system, where the utility is one of the four system utility functions introduced in Section |IIj 



Theorems 3.3 and 3.4 establish the strong NP-hardness of the problem for the special case c = 2 
and the general case c > 1, respectively. 

Theorem 3.3: The system utility maximization problem ([6]) with H = Hi, H2, H^, or H^ is 
strongly NP-hard when N/K = 2. 

Proof: Assume that N/K = 2. For any instance of the 3-dimensional matching problem, 



we construct the same system as in the proof of Theorem 3.1 and set pk = 5 {k E IC 
Theorem 3A\ directly implies that the following problem 



max min|i?fc| 

Pfc > Pfc > 0, keJC, neX, 

PkP] = 0' ^ J T^k, neU 
is strongly NP-hard, since the problem of checking whether its optimal value is greater than or 
equal to 3 log 2 is strongly NP-hard. 

Next, we prove that the sum-rate maximization problem 

1 ^ 

s-t. Y.n^MPl<Pk, keIC, (21) 
Pfc > Pfc > 0, ke)C, neAf, 
pIp^ = 0, W j j^k, neAf 
is also strongly NP-hard. In particular, we show that checking the optimal value of problem 



(21 1 is greater than or equal to 3 log 2 is strongly NP-hard. To this aim, consider the following 
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relaxation of problem pT]) 



max > R'u 

s-t- Efce/c EneATPfc < EfcGx; Pk = ^K, (22) 

Pfc > 0, kelC, neAf, 

PkP] = 0, V J ^ ^, neAf, 



where 



Due to the OFDMA constraint, we see that problem (22) is equivalent to 



max 



> 0, n G AT. 



(23) 



Noticing that problem ( [23] ) is convex, we can obtain its optimal solution 



^ (3, 3,. ..,3, 2, 2,. ..,2) 

and its optimal value 37^ log 2. Therefore, the optimal value Hi (i?*^ , R*2^ , R*^^ ) of the original 
problem ( [21] ) is less than or equal to 3 log 2, and the equality holds if and only if 

i?l^=i?;^ = ... = i?;,^=31og2. 

The latter holds if and only if the answer of the 3-dimensional matching problem is yes. 
Therefore, the optimal value of problem ( [2T] ) is greater than or equal to 3 log 2 if and only 
if the answer of the 3-dimensional matching problem is yes. 

For the cases H2 and H^,, notice that for all R2, . . . , Rk > 0, 

Hi{Ri, i?2, ■■■,Rk) < Hs{Ri, i?2, Rk) < -f^2(-Ri, R2, Rk) < Hi{Ri, R2, ...,Rk) 

and the equalities hold if and only if Ri = R2 = ... = Rk- Therefore, the optimal value of 
H2iRl^,R*2^, ...,R*kJ or H^{R\^,Rl^, ...,i?^J is greater than or equal to 3 log 2 if and only if 
the answer to the 3-dimensional matching problem is yes. This implies the strong NP-hardness 
of all the four utility maximization problems. ■ 
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Now consider the general case that c > 1. By constructing some dummy users and subcarriers 
as in the proof of Theorem 3.2[ we can extend Theorem 3.3 to the case c > 1 (Theorem 3.4). 
Take the sum-rate maximization problem as an example. It is simple to check that 

- in case 1 < c < 2, the sum-rate utility function of the constructed system is greater than or 
equal to 

3(c- l)log2 + (2 -c) log 11 

if and only if the given instance of the 3-dimensional matching problem with size (c— 1)K 
has a positive answer; 

- in case c > 1, the sum-rate utility function of the constructed system is greater than or 
equal to 

3{K - 1) log 2 + ((c - 2)K + 2) log 11 
K 

if and only if the given instance of the 3-dimensional matching problem with size K — 1 
has a positive answer. 

Hence, the sum-rate maximization problem in the case that c > 1 is strongly NP-hard. Similar 
results hold true for the other utility functions. 

Theorem 3.4: Given any constant c > 1, the system utility maximization problem Q with 
H = Hi,H2,H3, or if 4 is strongly NP-hard when N/K = c. 

We remark that the results in Theorems |3.3| and |3.4| are different from the one in It was 
shown in [9] that the sum-rate utility maximization problem (|6]) (i.e., H = Hi) is NP-hard when 
the number of users is equal to 2. The proof is based on a polynomial time transformation 
from the equipartition problem [16|, which is known to be NP-complete but not strongly NP- 



complete. Theorems |3.3| and |3.4| show that the sum-rate utility maximization problem is strongly 
NP-hard as the number of users increases. The other three utility functions (H = H2, H^, or 
H4) are not considered in [[9|. The reference [Q proved the NP-hardness of the three utility 
maximization problems (H = H2, H^, or H4) in the two-user case by establishing a polynomial 
time transformation from the equipartition problem. However, the proof of Q is nonrigorous, 
since the given equipartition problem has a yes answer does not imply that the transmission rate 
of the two users (at the solution of the corresponding utility maximization problem) is equal 
to each other, and vice verse. The complexity status of problem ([6]) with proportional fairness 
utility, harmonic mean utility, and min-rate utility remains unknown when K > 2 is fixed. 
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IV. Easy Cases 

In this section, we identify some easy cases when problem ([3]) or problem (|6]) can be solved 
in polynomial time. Before doing this, we introduce a concept called strong polynomial time 
algorithm. A problem is said to admit a strong polynomial time algorithm if there exists an 
algorithm satisfying the following two conditions: 

- the complexity of the algorithm (when applied to solve the problem) depends only on the 
dimension of the problem and is a polynomial function of the dimension; 

- the algorithm solves the problem to global optimality (not just e-global optimality). 

We remark that so far we do not know whether there exists a strong polynomial time algorithm 
to solve the general linear programming. When the interior-point algorithm is applied to solve 
the linear programming, it is only guaranteed to return an e-optimal solution in polynomial time 



and the complexity of the interior-point algorithm depends on the factor log(l/e) [22|. The best 



complexity results of solving the linear programming is still related to the condition number of 



the constraint matrix [23 1. 



A. Polynomial Time Solvability of Problem ([3]) when K = 1 

When there is only one user (i.e., K = 1), problem (|3]) becomes 



p = min > 



S.t. E„eArlog(l + ^) >7 



^Wo, (24) 
> > 0, n G TV. 



We can claim that solving problem ([24]) is equivalent to finding a minimal p such that the 
objective value of problem 



a p' 

max " ■ ■ ' 



^log 1 + 



S-t- EneArP"<P, ^^^^ 

> > 0, neJ\f 

is equal to 7. This is an important observation towards obtaining the closed-form solution of 



problem (24 1. In fact, if the optimal value of problem (24) is p, then the optimal value of problem 



( |25| ) is 7; and if the optimal value of problem ( [25] ) is 7, then the optimal value of problem ( [24] 
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is p. In more details, for any fixed p > 0, problem ( [25] ) is strictly convex with respect to {p"} 
and hence has a unique solution. Therefore, the set defined by 



1 + 



> 



7, 



(26) 



> > 0, n G 



contains only one point, which must be the solution of problem ( [25] ). Hence problems ( [24| ) and 



(25) share the same solution. 



The solution to problem ( [24) ) or problem ([25]) is given by the following extended water-filling 
solution 



P" 



mm < p 



1]" 



n 



eAf, 



(27) 



where A is chosen such that the objective value of problem ([25]) z^' equal to 7, and {x) 



max {x, 0} . After obtaining the optimal A, the optimal value p of problem ( [24] ) is given by 

1 7]"^ 

A 



mm < p", 



+ 



For completeness, we show in Appendix [A] that ( [27[ ) is indeed the solution to problem ([25]) and 
A is actually the Lagrangian multiplier corresponding to the constraint < p. 



We point out that the water-fiUing solution ([27]) extends the conventional water-filling solution 

1 r]" 
A a' 



P" 



mm 



n 



eAf 



(28) 



in [21 1 in the following two respects. First, the conventional water- filling solution ( [28] ) solves 
the power control problem ([25]) without the power budget constraints per subcarrier, while the 



power control problem ( 125] ) not only involves the total power constraint but also involves the 
power budget constraints per subcarrier. Second, the parameter A in the conventional water-filling 
solution ([28]) is chosen such that ^„g^p" = P, while the parameter A in ( |27] ) is chosen such 
that 



neAf 



1 + 



7- 



The only left problem now is to find r := 1/A in ( |27] ) such that the objective value of problem 
^5] ) is equal to 7. A natural way to find the desired r* is to perform a binary search on r, since 



the objective function of (25) is an increasing function with respect to r. As is known, the 
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efficiency of tlie binary search depends on the initial search interval of r. To derive a good 
lower and upper bound, we first order the sequence 



. (29) 



a" a' 



n=l 



and without loss of generality, suppose that 

Notice that it takes 0{2N log(2A^)) operations to order {6"}^^^^ . Then we calculate the objective 
values of problem ( [25] ) at r = 6" (n = 1,2, ...,2A^), denoting by {v'""}^^^. It follows from the 
monotonicity that 

<v^ <...<v''''. 

If there exists an index n such that t>" = 7, then r* = 6". Otherwise, we have t>" < 7 < t>"+^ for 
some n and hence r* G (6", 6"+^). In this case, we can start the binary search from (6", 6"+^) 
and it takes at most log ((6"+^ — 6")/e) iterations to obtain an e-optimal r. 

B. Strong Polynomial Time Solvability of Problem ([6]) when K = 1 

If the system has only one user, all the four system utility functions coincide and problem ([6]) 
becomes 

max > log 1 H 

{p"} ^ \ r ) 

s-t. E„,eArP"<P, ^^^^ 
p'^ > > 0, n G A/". 



The main difference between problem (31 1 and problem ([25]) lies in that, p is a given constant 



in problem (31), but is an unknown parameter in problem (25). This feature of p in problem 



ST]) enables us to design a strong polynomial time algorithm for the problem. 



Notice that the solution to problem ([3T[) is given by the extended water-filling solution ( [27] ), 
where r := 1/A is chosen such that = P- The desired f here can be found by solving 

a univariate linear equation. In fact, in a similar fashion, we can order {6"}^^^ in ([29]) and 



assume without loss of generality that ( [30[ ) holds. Then we calculate the total transmission power 
SneA/'^'" at T = {n = 1,2, ...,2A^), denoting by {t*"}^^!- It follows from the monotonicity 
that 
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If there is an index n such that p = u", then r = 6". Otherwise, we have that vJ^ < p < 
for some 1 < n < A^. Then for each n = 1, 2, A^, we have 



mm <; p ' I ^ " ^ 



p", if ^ + < 6*^; 



a' 



a" a' 
0, if — > 6". 



if — < 6" < — + 



(32) 



a" 



Therefore, the problem of finding the desired f reduces to solve a univariate linear equation in 
terms of r, and the desired f is obtained in a closed form. From the above discussion, we see 
that the complexity of finding the desired f is 0(2A^log(2A^)). 

We remark that Q has addressed problem pT] ), but without power budget constraints per 
subcarrier. 



C. Strong Polynomial Time Solvability of Problem ([3]) when N = K 

When the number of users is equal to the number of subcarriers (i.e., = K), we transform 
problem ([3]) into an assignment problem for a complete bipartite graph with 2K nodes in 
polynomial time. Specifically, we construct the bipartite graph Q = (V, E) with 

• the node set 

V^ = {1,2,...,K, l',2',...,K'}, 

where {1, 2, K} and {!', 2', K'} correspond to the set of users and the set of subcar- 
riers, respectively; 

• the edge set 

E = {{k,n) I A;G{l,2,...,ir}, n e {l\2' , K'}} 



with the weight w'l of edge {k,n) 

(exp (7fc) - 1) ril 



a 



if log ( 1 + 1 > Ik. 



k,k 



K K' 

X^X^^fc' Otherwise. 

k=l n=V 



(33) 
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Therefore, problem ([3]) can be equivalently reformulated as 



K K' 



{K) } = argmin 2^2^ 



w'Cx 



{_^k} k=l n=V 
K 



k-^k 



s-t. Eti< = l> n = V,2',...,K', (34) 
J2n=i' ^k — k = 1,2, K, 
xle{0,l}, k = l,2,...,K, n=V,2',...,K'. 
In the above problem, the binary variable is equal to 1 if user k transmits power on subcarrier 
n and otherwise. The first constraint Yl!k=i x^k — ^ — 1'' 2', K') stands for the OFDMA 
constraint, which requires that at most one user is allowed to transmit power on each subcarrier. 
The second constraint ^^=1' -"^k = 1 = 1)2, ■■■,K) requires that each user must transmit on 
one subcarrier to satisfy its specified transmission rate demand. 



From [20. Theorem 11.1], we know that the Hungarian method solves problem ( [34| ) to global 
optimality in 0{K^) operations. The Hungarian method is in essence a primal-dual simplex 
algorithm for solving the linear program 

K K' 

{{xD*} = arg min ^ ^ 



{^k} k=l n=V 



s-t. Eti< = l, n = V,2',...,K', (35) 
En=i' ^fe = 1' k = 1,2, K, 
Q<xl<l, k = l,2,...,K, n = l',2',...,K', 



which is a relaxation of problem ([34]) by replacing x^ G {0, 1} with < < 1. In general. 



the primal-dual simplex method is not a polynomial time algorithm. When it is used to solve 



problem ( [35] ), however, it can return the global optimal integer solution of problem ( |35] ) in 
0{K^) operations pol. 



We remark that if the optimal value p* of problem (34) satisfies 



K K' 
k=l n=V 



the optimal solution to (|3]) is 



[ply = {xlr^^^^^^^^^, k = l,2,...,K, n = l',2',...,K'- 
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Otherwise, we will have that 

K K' 
fc=l n=V 

which indicates that problem (|3]) is infeasible. 

In fact, problem Q is also polynomial time solvable when N = K + C, where C > 1 is a 
given constant integer. Specifically, we can first partition K + C subcarriers into K nonempty 
subsets {N'n]n=v ■ Denote the number of ways to partition K + C subcarriers into K nonempty 
subsets by S{K + C,K). We show in Appendix |b] that S{K + C,K) is upper bounded by 
{K + . After the partition of the subcarriers, the problem reduces to the case N = K. The 
only difference here is that the parameter (A; = 1, 2, n = 1', 2', -f^') is set to be 



= min 



s.t. 2^ log 1 + 1 > 



\ ^ (36) 



Vk 



Pi>pi>0, leAfn 
if the above problem is feasible; otherwise it is set to be 



K N 
k=l 1=1 



We know from Subsection IV-A that problem ([36]) is polynomial time solvable. Moreover, if 
problem ([36]) is feasible and J\fn = {/} only contains a singleton, we have 

(exp {-fk)-l)vi 



Then, for a given partition j9 = 1, 2, S{K + K), we can use the Hungarian method to solve 
problem ([34]) in 0{K^) operations, and assume the optimal value to be fp. Therefore, the optimal 



value of the original problem is 



min \v„} 

p=l,2,...,S{K+C,K) 



if it is strictly less than J2k=i J2^=iPk'^ otherwise the original problem is infeasible. It follows 
from the above analysis that problem ([3]) with = K+C can be solved in 0{K'^'^'^^) operations. 



December 21, 2012 



DRAFT 



25 



D. Polynomial Time Solvability of Sum-Rate Maximization Problem ([6]) without Total Power 
Constraint 

If we drop the total power constraint XlneA/'^'fc — Pk (k ^ K,), then problem ([6]) with H = Hi 
becomes 

max Hi(Ri,R2,...,Rk) 

s.t. > > 0, keJC, n G AT, (37) 
pIpI = 0, y 3 7^k, neU. 



Problem ( [37| ) is always feasible as = (for all k and n) is a certificate for the feasibility. 
Further, since problem ( [37] ) does not involve the total power constraint for each user, its optimal 
solution (k E IC, n E M) is either or p^. To solve problem ([37]) in polynomial time, we 



consider to transform it into the polynomial time solvable Hitchcock problem [20| (also known 
as transportation problem). 

Hitchcock problem. Suppose there are K sources of some commodity, each with a supply 
of Ofc E (k = 1,2, ...,K) units, and terminals, each of which has a demand of d„ E 
]R+ (n = 1, 2, ...A^) units. Suppose that the unit cost of transporting the commodity from source 
k to terminal n is E (k = 1,2, K, n = 1, 2, ...,N). The problem is how to satisfy the 
demands at a minimal cost? 

In fact, by setting 

ai, = N ik = 1,2,..., K); 
dn = l{n = l, 2, ...N), dN+i = {K- 1)N; 

( Oik kPk \ 

4 = c-\og\^l + ^^j >0, k = l,2,...,K, n = l,2,...,N, 
c^+i = 0, k = l,2,...,K, 

^'LkPk 



where c = maxfc_„ <^ log ( 1 H — ) |> , we can see that problem ( [37] ) is equivalent to the 
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following Hitchcock problem 

K N+1 
l^fcj k=l n=l 

s-t- Ef=i^fc = 1, r2 = l,2,...,iV, 

ELa:f+^ = (^-W ^^^^ 
En=l^^ = A^, fc = l,2,...,if, 
a:^ = {0,l}, A; = l,2,...,ir, n = 1, 2, iV. 
In the above problem, the binary variable = 1 if user k transmits full power on subcarrier n 
and = if user k does not transmit any power on subcarrier n. The constraint Ylk=i ^fc ~ -'- 
and the constraint x1 = {0,1} (k = 1,2, ...,K) implies that at most one user is allowed to 
transmit power on each subcarrier n. However, one user is allowed to transmit on multiple 
subcarriers. The variables x^~^^ {k = 1,2,...,K) are dummy variables, since c^^^ = for all 
k = 1,2,...,K. Moreover, we know from [20, Theorem 13.3 and its Corollary] that problem 



( [38] ) is equivalent to the hnear program 

K N+1 

{(4r} = min J^J^cX 



{^k] k=l n=l 

s-t- Ef=ia:^ = l, n = l,2,...,N, 

T.n=i^l = N, k = l,2,...,K, 
0<a;^<l, k = l,2,...,K, n = l,2,...,N. 



(39) 



The equivalence between problems (38) and (39) is because the linear equality constraint in (39) 



satisfies the so-called totally unimodular property [20|, and thus all the vertices of the feasible 



set of problem ( [39] ) are integer. Since the linear program ( [39[ ) is polynomial time solvable, 
the sum-rate maximization problem ([37]) is polynomial time solvable. Further, if {(x^)*} is the 
optimal solution to problem (|39]), the optimal solution {(p^)*} to problem ( fTTj ) is 

{pir = KrPl. k = l,2,...,K, n = l,2,...,N. 
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In a similar fashion, we can show that the weighted sum-rate maximization problem 

K 



max WkRk 

(40) 



^^^^ fc=i 



s.t. > > 0, e /C, n G AT, 
PkP] = 0, V j ^ A;, n e 
is also polynomial time solvable, where Wk > {k = 1,2,...,K) are non-negative weights 



Again, we can transform problem (40) into the Hitchcock problem in polynomial time. The 



corresponding parameters of the Hitchcock problem are the same as before, except 

cl = c-Wk\og[l + ^^j >0, k = l,2,...,K, n = l,2,...,iV, 

'^k,kPk 
Vk 



where c = max^ „ < Wk log 1 + 



We use the following theorem to conclude Section IV 
Theorem 4.1: The following statements are true. 

- Problem (|3]) is strongly polynomial time solvable when the number of users is equal to the 
number of subcarriers. 

- Problem ([3]) is polynomial time solvable when there is only one user in the system. 

- Problem ([6]) is strongly polynomial time solvable when there is only one user in the system. 

- The (weighted) sum-rate maximization problem ([6]) without the total power constraint is 
polynomial time solvable. 

V. Concluding Remarks 

Dynamic allocation of subcarrier and power resources in accordance with channel and traffic 
load changes can significantly improve the network throughput and spectral efficiency of the 
multi-user multi-carrier communication system where a number of users share some common 
discrete subcarriers. A major challenge associated with joint subcarrier and power allocation is 
to find, for a given channel state, the globally optimal allocation strategy to minimize the total 
transmission power or maximize the system utility function. This paper mainly considers the 
computational challenges of the joint subcarrier and power allocation problem for the multi-user 
OFDMA system. We have shown that the general joint subcarrier and power allocation problem 
for the multi-user OFDMA system is strongly NP-hard as the number of the users in the system 
increases. The complexity result suggests that we should abandon effort to find globally optimal 
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subcarrier and power allocation strategy for the general OFDMA system unless for some special 
cases (i.e., the case when there is only one user in the system, and the case when the number 
of users is equal to the number of subcarriers). The problem is shown to be polynomial time 
solvable or strongly polynomial time solvable in those special cases. In a companion paper, we 
shall design efficient algorithms to solve the joint subcarrier and power allocation problem. 
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Appendix A 



The KKT condition of problem ( [25] ) is 
( -1 



- + A + ^" - z/*^ = 0, 



(41) 



E„,eArP" < A > 0, A(E„eArP" - P) = 0, 

> 0, z/" > 0, p"z/" = 0, n G TV, 
where A, i/" are the Lagrangian multipliers corresponding to the constraints XlnGA^^*" — ^' 



> p", and > 0, respectively. Recall that problem ( [25| ) is convex for any fixed p > 0. It 
follows that, to show ([27]) is the solution to problem ([25]), it suffices to show that there exist 



appropriate Lagrangian multipliers such that (p", A, i/") satisfy the KKT system ( [4T] ). 

Specifically, we choose A > such that the objective value of problem ( |25| ) is equal to 7, and 



if = in ([27]), then 1/A - ifjd^ < 0; we set ^" = and z/" = A - a"/r7" > 0; 
if < p" < in ([27]), then > 1/A - r^^/a" > 0; we set ^" = z/" = 0; 
if = in ([27]), then 1/A - r/"/a'' > p"; set C = — — - \>0 mdu" 



0. 



It is simple to check (p"^, A, z/") constructed in the above satisfies the KKT system ( [4T] ). Hence, 
( |27| ) is the solution to problem ( |25| ). 
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Appendix B 



In combinatorics, a Stirling number of the second kind [24|, denoted by S{N,K), is the 
number of ways to partition a set of objects into K nonempty subsets, where > i^' > 1. 
Stirling numbers of the second kind obey the following recursive relation 

S{N + 1,K) = S{N,K -1) + KS{N,K) (42) 

with initial conditions S{k, k) = 1 and S{k,0) = 1 for any > 1. To understand this formula, 
observe that a partition of 1 objects into K nonempty subsets either contains the {N + l)-th 



object as a singleton (which corresponds to the term ^(A^, K — 1) in (|42])) or contains it with 
some other elements (which corresponds to the term KS{N, K) in (|42])). 

Next, we claim by induction that S{K + C,K) < {K + C^^ . In fact, when C = 1, we can 
obtain 

This is because that dividing K + 1 elements into K sets means dividing it into one set of size 
2 and K — 1 sets of size 1. Therefore, we only need to pick those two elements. Assume that 

S{K + C ~l,K) < {K + C - lf'^^'^\ we show that S{K + C,K) < {K + Cf^. Invoking 
( |42l ), we have 

S{K^C,K) = S{K + C -1,K -1) + KS{N + C -1,K) 

< S{K + C-1,K-1) + {K + Cf^-^ (43) 

< {K + cy^, 

which shows that S{K + C,K)<{K + Cf^ holds true. 
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